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The accuracy of a single s-orbital representation of Cu towards enabling multi¬ 
thousand atom ah initio calculations of electronic structure is evaluated in this work. 
If an electrostatic compensation charge of 0.3 electrons per atom is used in this basis 
representation of copper, the electronic transmission in bulk and nanocrystalline Cu 
compares accurately to that obtained with a Double Zeta Polarized basis set. The use 
of this representation is analogous to the use of single band effective mass represen¬ 
tation for semiconductor electronic structure. With a basis of just one s-orbital per 
Cu atom, the representation is extremely computationally efficient and can be used 
to provide much needed ab initio insight into electronic transport in nanocrystalline 
Cu interconnects at realistic dimensions. 
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Copper is the current and and projected metal of choice for back-end-of-line (BEOL) 
interconnects used in semiconductor logic and memory technology. It has been well known 
that resistivity trend for copper is strongly non-linear and has been increasing rapidly with a 
reduction in interconnect width. A number of semi-empirical models have been formulated 
to explain this phenomenon. Broadly speaking, these models can be categorized as variants 
of the Fuchs-Sondheimer model (F-S)P21of surface scattering-induced increase in resistivity 
or the Mayadas-Shatzkes (M-S) modeP of grain-boundary scattering-induced increase in 
resistivity. 

While these models provide important insights into electron transport in metals, they are 
£t a posteriori to experimental data with empirical parameters such as average reflectivity 
and specularity. Such parameteric models can be £t with a non-unique set of parameters 
and consequently provide very little predictive insight from a materials design perspective. 
For instance, new experimental data on nanocrystalline copper requires a recalibration of 
the model to account for the new data. Often times, such empirical models can lead to 
conflicting physical insight. For instance, Graham et al.^ht the resistivity of sub-100 nm Cu 
interconnect lines to a purely surface scattering based model, while Steinhogl et al.^ £t the 
resisitivity of Cu lines in the same dimensional range to a predominantly grain boundary 
scattering based model. 

In spite of the apparent seriousness of the resistivity runaway problem and lack of un¬ 
derstanding of the fundamental features that govern electronic transport in nanocrystalline 
(polycrystalline, with nanometer-range sized grains) Cu structures, there have been very few 
hrst principles based investigations on electron transport in nano crystalline Cii^®®. Part of 
the reason for such a paucity of hrst principles based models is the cumbersome computa¬ 
tional requirement (memory and execution time) for convergence of hrst principles transport 
calculations for realistic Cu nanocrystalline systems. This is illustrated with a simple exam¬ 
ple - The ITRS projects the smallest damascene interconnect cross sectional area of 10 x 22 
nm^ for damascene Cu in the year 202^^. Cu has a lattice constant of 3.61 A at room tem¬ 
perature. A cross sectional sliver of a length of just 1 nm with the aforementioned area of 
Cu at the projected ITRS dimensions contains approximately (assuming bulk density holds 
even for nanocrystalline samples) 18700 atoms. If one were to simulate such a sliver using 
linear scaling order-N Linear Combination of Atomic Orbitals (LCAO) density functional 
theory (DFT) calculations with a moderately sized basis set (Double Zeta Polarized (DZP)), 
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containing 2 s, 10 d and 3 p orbitals, one would require approximately 1.25 TB to store a 
full representation of this system in memory (64 bit word size) for DFT calculations. Even 
if one were to reduce this to a sparse representation, one would require several hundred 
GB since s orbital interactions in Cu are long ranged. Apart from the fact that memory 
requirements (for full storage) for a system of size N atoms and an M orbital basis per 
atom scale as {N x M)^, full eigenvalue solves (required for density computation in metals) 
scale as 0((A^ x M)^). This is in addition to the 0(A^ x M) time complexity of evaluating 
Overlap and Hamiltonian integrals in DFT)42!. 

What is needed then, is an ab initio method that is computationally efficient for systems 
containing several thousand Cu atoms while retaining an acceptable level of accuracy. In 
this respect, Cu possesses a unique advantage as compared to other metals such as partially 
d hlled transition metals. Cu has a valence electron conhguration of The d orbitals 

are spatially localized and completely filled while the s orbital is only partially filled and 
delocalized. An obvious question then arises - is it feasible to to represent Cu’s electronic 
structure ab initio using only the partially hlled and delocalized s orbitals while retaining 
acceptable accuracy? If the answer were to be in the affirmative, such a method would allow 
tremendous computational savings, allowing the ab initio investigation of multi-thousand 
atom Cu systems. 

In this paper we report the results of our investigation into this question. We provide 
conditions for which the electronic structure of systems represented in a single s orbital 
representation can be compared favorably with that of DZP representations, which have in 
turn compared favorably to experiment in previous workP^^^HI g^d goal in this investi¬ 
gation is evaluation of the feasibility of using a reduced basis in studying electron transport 
in nano crystalline Cu structures where the scattering mean free path is signihcantly smaller 
than the phonon scattering mean free path. Accordingly, we use ballistic transmission and 
ballistic transmission per unit area as figures of merit for comparing electronic structure 
throughout this paper. 

The rest of this paper is organized as follows. We first describe the computational method¬ 
ology used in our investigation. This is followed by a comparison of electronic structure 
computed using the two sets of bases for a variety of boundary conditions and polycrys¬ 
talline configurations. We then conclude with a discussion on the relative accuracy of the 
method an its applicability in transport calculations on multi-thousand atom systems. 
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Double Zeta Polarized basis set 


Single 4s orbital basis set 




FIG. 1. Radial part of the Numerical Orbitals used in the DFT calculations. On the left is a 
Double Zeta Polarized (DZP) basis set containing s, p and d orbitals. On the right is a single 
orbital basis set containing just one s orbital. 

We performed first principles LCAO DFT calculations using the Atomistix Toolkit Pack¬ 
age (ATK)II^ within the Local Density Approximation (LDA) and the exchange correlation 
functional of Perdew And Zunger (PZ)II^. The reference basis was a numerical DZP basis 
consisting of s orbitals, d orbitals and p orbitals. The radial part of this basis is shown in 
hgure [Tj The radial part of the single 4s orbital numerical basis is also shown on the right 
in the same figure. 

The band structure and density of states (DOS) computed using the two bases and bulk 
k grid of 20^ k points is shown in figure An important feature of Cu’s electronic structure 
is the ’neck’ at the L symmetry point in the Brillouin Zone. This neck is situated about 
1.23 eV below the Fermi Energy in the DZP computed band structure. If the s orbital basis 
is used as is, the neck occurs 1.3 eV above the Fermi Energy. The bands must be aligned 
correctly to reproduce key features of the Fermi Surface with acceptable accuracy. In order 
to do this, we compensated (electrostatically ’doped’) the bulk Cu configuration with 0.3 
electrons. This is equivalent to introducing a overall negative compensation charge of 0.3 
electrons per atom in the unit cell. The resultant band structure reproduces the relative 
position of the L neck as shown in hgure 

It is also evident from hgure that the curvature of the bands in the s orbital basis is 
signihcantly higher than that computed in the DZP orbital basis. On the other hand the 
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Bandstructure Density of States 




FIG. 2. Comparison of band structure (left) and DOS (right) computed in LDA DFT using the 
two bases. 


DOS in the s orbital basis is significantly lower. Given that the net ballistic transmission 
is a product of DOS and average group velocit}^, it can be expected that the increase in 
group velocity is balanced out by the increase in DOS and that the ballistic transmission 
computed using the two bases is comparable. This is evident in the k-resolved transmission 
spectra computed in the Non Equilibrium Greens Function (NEGF) formalismP^n^ shown 
in hgure 
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FIG. 3. Transverse k resolved Transmission spectra for Cu computed using the DZP basis (top) 
and single s-orbital basis (bottom) 

The transmission spectra shown in hgure were computed at the Fermi Level for the 
smallest unit cells in [100], [110] and [111] transport orientations. While it is evident from 
this hgure that while the Fermi transmission spectra for the diherent transport orientations 
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computed using the single s-orbital basis is not exactly the same as that computed using 
the DZP basis, the qualitative similarity between the two is quite signihcant considering 
that the s-orbital basis set is a severely reduced representation. In order to effect a more 
quantitative comparison, the transverse momentum-averaged transmission was computed 
for each energy point and a net transmission was obtained by Fermi-averaging the energy 
resolved transmission. This was then done for multiple cross sectional areas. The resultant 
Transmission versus Cross Sectional Area plot is shown in figure 


Bulk Nanowire 



FIG. 4. Transmission versus Cross Sectional Area for bulk and nanowire copper plotted for different 
transport orientations Computed using the s-orbital and DZP bases 

It is evident from the computed bulk transmission and Fermi Surfaces that the s-orbital 
model matches the DZP results quantitatively in bulk. More interesting for practical ap¬ 
plications, however, is the behavior under confinement and when sources of disorder such 
as grain boundaries are present. To facilitate this comparison, Cu nanowire cross sections 
oriented along different orientations with a 1 x 2 aspect ratio were constructed. The average 
transmission along transport direction was computed versus cross sectional area and plotted 
alongside the bulk values in figure It is evident that the nanowire transmission computed 
using the two sets of bases matches well quantitatively. The trend in transmission mismatch 
between orientations for different orientations is also matched. 

To compare the transmission in grain boundary structures, structures oriented along 
the canonical orientations [100], [110], [111] were combined and four separate boundary 
conditions were applied 
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Bulk (periodic in transverse direction) with a single grain boundary and open (device) 
boundary conditions along transport direction. 


• Bulk with grain boundary repeated inhnitely along transport direction. 

• Nanowire with a single grain boundary and open boundary conditions along transport 
direction. 

• Nanowire with grain boundary repeated inhnitely along transport direction. 

Since the number of possible conhgurations considering grain size and orientation distribu¬ 
tion is potentially inhnite and computationally cumbersome to compute in the DZP basis, 
we limit ourselves to a limited number of grain orientation distributions and a grain size of 
Inm. 


[100]/[110] Bulk Single Grain Boundary ...[111]/[111]... Bulk Infinite Twin Boundaries ( 13 ) 





FIG. 5. Examples of the grain boundary structures for which transmission was computed in the 
s-orbital and DZP bases. The black lines indicate the unit cell used in the calculation. For single 
grain boundary structures the area on either side of the grain boundary is repeated semi-infinitely 
in the transport direction. 


An example of a structure for each boundary condition is shown in hgure Tables [I] and 
[II] compare the transmission per unit area computed using the two sets of bases. Of special 
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interest is the S3 [111]/[111] twin boundary that has the lowest specific contact resistivity 
of all grain boundaries and is the least disordered grain boundary structure in Cu. 

TABLE I. Comparison of Single Grain Boundary Transmission per Cross Sectional Area for Double 
Zeta Polarized (DZP) and s-orbital basis 


Bulk Nanowire 


Grain Distribution 

T/A (A-2) 

Grain Distribution 

T/A (A-2) 


DZP 

s-orbital 


DZP 

s-orbital 

[100]/[110] 

0.119 

0.113 

[100]/[110] 

0.092 

0.096 

[100]/[111] 

0.118 

0.112 

[100]/[111] 

0.075 

0.074 

[110]/[111] 

0.112 

0.101 

[110]/[111] 

0.063 

0.070 

[111]/[111] E3 

0.135 

0.145 

[111]/[111] S3 

0.102 

0.099 


The results in figures through and tables through show unambiguously that 
the transmission computed using the single s-orbital basis matches the qualitative trends 
in the transmission computed using the DZP basis, regardless of boundary condition, grain 
orientation distribution or structure. Quantitatively speaking, an excellent match is obtained 
for bulk transmission, monocrystalline and grain boundary structures. For nanowires, an 
extremely accurate match is obtained for monocrystalline structures. For grain boundary 
structures, however, the quantitative match is not as accurate even if the qualitative trend is 
accurate. For instance a maximum error of 36 % is obtained for the the nanowire structures 
(multi grain boundary [100]/[110] case). 


TABLE If. Comparison of Multiple Grain Boundary Transmission per Cross Sectional Area for 
Double Zeta Polarized (DZP) and s-orbital basis 



Bulk 


Nanowire 


Grain Distribution 

T/A (A-2) 

Grain Distribution 

T/A (A-Q 


DZP 

s-orbital 


DZP 

s-orbital 

...[100]/[110]... 

0.079 

0.082 

...[111]/[100]... 

0.033 

0.021 

...[100]/[110]/[111]... 

0.070 

0.053 

...[100]/[111]... 

0.024 

0.019 

...[111]/[111]... S3 

0.130 

0.140 

...[111]/[111]... S3 

0.076 

0.083 










Nevertheless, there are important reasons why this mismatch should not prevent the use 
of the single s-orbital basis in transport calculations in lieu of DZP or other ’full’ basis sets. 
The cross sectional areas considered in this paper were limited to a maximum cross sectional 
area of 181 owing to difficulties in converging structures larger than 800-900 atoms using 
the DZP basis set. This is far smaller than current or projected interconnect dimensions. It 
is evident from the results above that as the surface area to volume ratio increases (as one 
goes towards bulk systems), the accuracy of transmission calculated in the s-orbital basis 
set improves. Additionally, it is important to note that the projected resistivity runaway 
trencP^ is non-linear and can be £t with exponential behavior. Even with an error bar of 
20 to 40 % and correct qualitative trends an immense amount of insight can be gained 
about resistivity behavior of Cu at current and future nodes. In addition, the s-orbital basis 
could be very useful in gaining relatively quick insight about the lower limits of Copper 
resistivity at a given dimension or the effects of surface roughness on Cu resistivity at larger 
cross sectional areas. An added advantage of using this basis is that there are no further 
adjustable parameters once the bulk bands are aligned using doping (in our case, a doping 
of 0.3 electron per atom was used throughout). 

It should also be pointed out that the feasibility evaluation carried out above is restricted 
to electronic structure. This is similar in spirit to attempts at modeling semiconductor 
electronic structure by use of a single effective mass. Total energy calculations using a 
single s-orbital give results that do not match DZP total energies since d orbitals contribute 
signihcantly to total energy in Cu. 

In conclusion, we have evaluated the accuracy of a single s-orbital representation of Cop¬ 
per towards enabling multi-thousand atoms ab initio calculations of electronic structure. 
We found that upon modihcation through doping, the s-orbital representation accurately 
reproduces key features of the bulk Fermi surface, the electronic structure of monocrystalline 
nanowires, bulk and nanowire single grain boundary structures and twin boundary struc¬ 
tures. The electronic structure of multi grain boundary structures showed less comparative 
quantitative accuracy, but accurately matched all qualitative trends obtained through DZP 
calculations. The use of this representation is analogous to the use of single band effective 
mass representation for semiconductor electronic structure. We have found that the repre¬ 
sentation allows easy scaling to multi-thousand atom systems due to use of a reduced basis. 
With a basis of just one s-orbital per Cu atom, the representation is extremely computa- 
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tionally efficient and can be used to provide much needed ab initio insight into electronic 
transport in nanocrystalline Cu interconnects containing several thousand atoms. 

We thank Mark Rodder for a careful reading of the manuscript and for helpful suggestions 
made during the course of this work. 
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